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ABSTRACT 


The results of an Investigation of the influence of simulated 
turbulence on aircraft handling qualities is presented. Pilot opinion 
of the handling qualities of a light general aviation aircraft were 
evaluated in a motion-base simulator using a simulated turbulence 
environment. A realistic representation of turbulence disturbances 
is described in terms of mis intensity and scale length and their 
random variations with time. The time histories generated by the 
proposed turbulence models showed characteristics which appear to 
be more similar to real turbulence than the frequently-used Gaussian 
turbulence model. In addition, the proposed turbulence models can 
flexibly accommodate changes in atmospheric conditions and be easily 
implemented in flight simulator studies. 

Six turbulence time histories, including the conventional Gaussian 
model, were used in an IFR- tracking task. The realism of each of the 
turbulence models and the handling qualities of the simulated airplane 
were evaluated. Analysis of pilot opinions shows that at approximately 
the same ms intensities of turbulence, the handling quality ratings 
transit from the satisfactory level, for the simple Gaussian model, 
to an unacceptable level for more realistic and compositely structured 
turbulence models. 
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CHAPTER I 


INTRODUCTION 

Simulated time histories of aircraft motion in a turbulence 
environment are required in a variety of engineering applications, 
and their use appears to be increasing as more intricate and 
sophisticated design studies are attempted. As an example, the use 
of flight simulators for the study of airplane handling and ride 
quality has proven to be more valuable when disturbances in the form 
of artificially simulated turbulence are Introduced into the system. 
Several methods have been used to generate turbulence signals; each 
one aimed at realizing the actual atmosphere as closely as possible. 

A realistic representation of turbulence becomes especially important 
in the simulation of future aircraft with high sensitivity to 
turbulence, as even light to moderate turbulence may seriously degrade 
their controllability and ride quality. Low altitude atmospheric 
turbulence critically effects the evaluation of vehicle handling 
qualities, pilot work load, ride quality, and other design factors. 
Several empirical studies (1,2,3) have shown that low altitude clear 
air atmospheric turbulence is only locally isotropic, i,e. , isotropic 
over a finite range of scale lengths. The proposed gust model accounts 
for the anisotropy of typical low altitude clear air turbulence by 
randomly varying the rms velocities and scale length of the gust field. 
The scale lengths predicted by either the Von Karman or the Dryden 
models (4) are large compared to real atmospheric turbulence and hence 
the scale length distribution is modified to achieve compatibility. 
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With a suitable combination of scale length and intensity 
distribution, the proposed model will simulate various atmospheric 
conditions characterized by altitude, stability, and terrain. This 
new model is mechanized to be included In a flight simulator experiment 
in order to determine to what extent the pilots are sensitive to 
changes in atmospheric conditions and the realism of the model. The 
following chapters describe the proposed turbulence model and the 
flight simulator experiment In detail. 



CHAPTER II 


LITERATURE SURVEY 

In this chapter statistical properties of atmospheric turbulence 
are reviewed and presently-used simulation techniques are discussed. 

A review of basic definitions in probability and statistics is 
included in Appendix A, 

2 . 1 Properties of Atmospheric Turbulence 

Simulations of aircraft flying through atmospheric turbulence 
require a realistic model of the physical environment. Therefore, 
simulation studies in general begin with a study of the real 
atmosphere. In references 1-3, atmospheric data have been reported 
characterizing various atmospheric conditions for variation in 
terrain, stability, altitude, temperature, time, season, and 
geographic location. This data has been suitably modified to 
establish a basis of comparison for the simulated turbulence field. 

The following criteria are used as the bases of comparison: 

a) Output Statistics 

Mean and standard deviations of the gust velocities. 

b) Probability Distribution 
Cumulative probability 
Probability density 

Fourth and sixth normalized moments 

c) Patchiness of the Field 

d) Power Spectral Density 

e) Element of Surprise 
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Each of these properties will be discussed from the standpoint of real 
atmospheric turbulence. 
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Mean ; Analysis of several sets of data presented in Reference (1) 
indicates that the mean velocity of atjnospheric turbulence is 0.0 + 0.1 
ft/sec (0 + 0.03 m/sec) . 

Standard Deviation ; The standard deviation of the velocity field 
for low altitude clear air turbulence is 3.0 + 1.31 ft/sec (0.91 + 0.4 
m/sec). Typical values for various conditions are listed in Reference 
(4) as : 

= 2 ft/sec (0.61 m/sec) for light turbulence 
0 ^ ■ 4 ft/sec (1.22 m/sec) for moderate turbulence 
0 ^ » 6 ft/sec (1.82 m/sec) for severe turbulence 
where u is the longitudinal gust component. 

Probability Distribution ; The probability distribution of a 
random process provides information concerning the range of values 
assumed by that function and the frequency with which they occur. 

As there is little experimental data available which distinguishes 
between probability distributions of different gust components, no 
distinction will be made here. 

Probability Density Distribution ;: Figure 1 presents data from 

Reference (5) showing a typical probability density distribution of 
atmospheric turbulence velocity. The departure from the Gaussian 
curve clearly Indicates Increased probabilities of large and small 
gusts. 

Normalized Central Moments ; The fourth and sixth moments of low 
altitude real atmospheric turbulence are = 3,5 and = 21.7, 


respectively. ( 5 ) 




Probability % 
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Gust Velocity Component 
w - ft/sec (m/sec) 


Figure 1 Typical Probability Density Function of Atmospheric 
Turbulence 
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Patchiness ; It is knovm that turbulence has a non-Gaussian 
patchy structure which seems to occur in bursts of relatively intense 
motion separated by areas of relative calm. Figure 2 shows typical 
patchy characteristics for a 40-sec sample of real atmospheric 
turbulence. 

Power Spectral Density (PSD) : The PSD of a random process 

provides information on the average contribution to the process from 

the frequency components which make it up. Figure 3 presents a 

typical plot of PSD of low altitude clear air atmospheric turbulence. 

It may be observed that at high frequencies, the spectral density 

-2 

varies as inverse square of frequency (m ). On the other hand, at 
low frequencies the PSD is characterized by a horizontal asymptote. 
Two convenient mathematical forms are used to represent the power 
spectra of atmospheric turbulence. These are; 
a) Von Karman Spectra 


CiJj) 
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( 2 . 1 . 1 ) 


( 2 . 1 . 2 ) 


(2.1.3) 


and 




Figure 2 Typical Patchy Mature of Atmospheric 
Turbulence (6) 
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(2.1. A) 


( 2 . 1 . 6 ) 


( 2 . 1 . 6 ) 


Uq = initial total velocity 

= scale for 1^^ turbulence velocity 

= nns gust intensities 
u = frequency 

1 = u, V, w gust components. 

The Von Karman spectral shapes, although accurate, are not con- 
venient for turbulence modelling work since they cannot be matched 
using linear filters- This is due to the noninteger power appearing 
In the denominators. Thus, in order to avoid computational complexity 
in this report, the Dryden form is adopted. 

Element of Surprise ; More often than not, real atmospheric 
turbulence, when encountered, presents an element of surprise. It is 
not easy to formulate a model of this phenomenon in terms applicable 
to flight simulator work. It seems that a measurement of "sudden jump" 
in the velocity field can be used as a possible criterion to describe 
this phenomenon. Relative frequency of "sudden jump" of atmospheric 
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turbulence can be compared to the simulated turbulence field. Changes in 
aircraft orientation angles can also be used to measure this phenomenon. 
2.2 Presently-Used Simulation Techniques 

In this section several presently-used simulation techniques are 
discussed from the standpoint of their statistical realism and suit- 
ability for use in flight simulators. 

Measured Turbulence Field ; Flight recordings of atmospheric 
turbulence is perhaps the most obvious method of producing a realistic 
simulation. There can be little argument as to whether or not these 
time histories are an accurate and realistic representation. However, 
it is difficult to adjust the measured time histories to allow for 
conditions other than those for which it was recorded. No allowances 
can be made for changes of altitude or different atmospheric conditions. 
Another serious drawback is that the recorded time histories are fixed 
in length. Extended run times, therefore, cannot be accommodated 
without repetition. From the simulation point of view, the pilots 
tend to recognize some of the characteristics of the turbulence field 
and develop an intuition for predicting the field. This defeats the 
purpose of an artificially simulated turbulence field, which is to 
provide unpredictable external disturbances. It can, therefore, be 
concluded that flight recordings of atmospheric turbulence are not 
suitable for the simulation of typical turbulence. 

Sum of Sine Waves : Reference (5) describes this method in summary 

form. This technique involves superimposing several sinusoidal waves 
of different frequencies and amplitudes. The resultant is used to 
represent time histories of turbulence. One obvious disadvantage of 
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this method is that it contains only a finite range of frequencies 
whereas actual atmospheric turbulence consists of an infinite number 
of frequency components. 

Results of this simulation are not available but the model can 
justifiably be discarded on the basis of its inadequacy in matching 
the frequency content. 

Method of Orthogonal Functions ; In this method (7) , the recorded 
time histories of turbulence are decomposed into eigenfunctions of a 
covariance matrix. The probabilistic structure of the eigenfunction, 
and the coefficients of each of the time histories are studied. 

Simulated time histories are then regenerated by suitably modifying 
the distribution of the coefficients. The available preliminary 
results show that this technique adequately models the frequency contents 
and also presents an element of surprise. However, this model fails to 
show a patchy non-Gausslan characteristic which is typical of the real 
atmosphere. In addition to the mathematical complexity of the 
technique, its application is limited since recorded time histories 
are needed. 

Gaussian Turbulence Model : The classical method, most widely used 

for turbulence simulation, is the linearly filtered white noise tech- 
nique. Here the turbulence gust field is produced by passing white 
noise through a linear filter as shown in Figure 4a. The resultant 
signal is shaped so that the power spectrum and rms intensities match 
those of real turbulence. A Dryden or Von Karman form (6) are normally 
used to model the power spectrum. This model is remarkably easy to 
Implement and can be adjusted for any general power spectrum. However, 
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Figure 4a 


Gaussian Turbulence Simulation 



40 Sec. Sample 


Gust 

Velocity 

Component 


Figure 4b Derivative of Vertical Velocity 
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this model too falls short of reproducing the non-Gaussian patchy 
nature of real turbulence. Figure 4b compares the artificially 
simulated gust field using the Gaussian model (with a Dryden spectrum) 
and real atmospheric turbulence. It may be observed that the Intensity 
for the Gaussian model is nearly constant whereas measured ("real") 
turbulence exhibits a patchy nature or Intensity bursts. Test pilots, 
when exposed to this model in a flight simulator, rated the realism 
fair to poor. (5,6) 

Non-Gaussian Turbulence Model ; Reference (6) presents a non- 
Gaussian turbulence model. Time histories are generated by multiplying 
two independent random variables, one to represent the turbulence 
within a patch and the other to represent the variation of intensity 
with time. Figure 5a shows two independent Gaussian white noise 
generators and linear filters, which produce Gaussian random variables, 
a(t) and b (t) . These variables are then multiplied to produce gust 
time histories. 

The non-Gaussian model proposed in Reference (5), a modification 
of the above, is shown in Figure 5b. Here a(t), b(t), and d(t) are 
independent Gaussian processes. The process c(t) is generated by 
multiplying a(t) and b(t) , The resultant process, c(t), a modified 
Bessel process. Is summed with d(t) to form the output, u(t). The 
most remarkable achievement of this model is that the patchy character- 
istic and several statistical parameters of the simulated turbulence 
field can be varied simultaneously by varying the standard deviation 
ratio (R ■ a^/a^) . However, when R is varied to achieve one set of 
statistical properties, several other statistical parameters of 
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Figure 5a Non-Gaussian Turbulence Model (Reference 6) 
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Figure 5b Non-Gaussian Turbulence Simulation (Reference 5} 








16 


interest do not match real turbulence. In addition, due to the 
mathematical complexity, the mechanization of this model on a flight 
simulator is complicated and expensive. 

It can be observed from the review of presently-used simulation 
techniques that there is a need for a new model which adequately 
matches real atmospheric turbulence and is simple to implement in 
flight simulator studies. None of the preceding models have the 
flexibility of simulating various atmospheric conditions characterized 
by altitude, stability, and terrain. It is, therefore, necessary to 
Introduce a new turbulence model which is realistic and can flexibly 
accommodate changes in atmospheric conditions and be easily Implemented 
in flight simulator studies. 



CHAPTER III 


PROPOSED GUST MODELS 


Of the simulation techniques described, the Gaussian turbulence 
model is the simplest to implement and least expensive computationally. 
The proposed turbulence models, modification of the Gaussian simulation 
technique, retain the simplicity of the Gaussian technique while 
adequately modelling the characteristics of real atmospheric 
turbulence. In this report three basic models are proposed: 

1) Modified Gaussian Model 

2) Rayleigh Model 

3) Variable Length and Intensity (VLI) Turbulence Model. 

3.1 Modified Gaussian Model 

A block diagram of the modified Gaussian model is presented in 
Figure 6. Gaussian white noise, <{>q, is passed through a linear filter, 
G^Cs) i “ u, V, w, whose power spectrum is given by a Dryden model 
(e.g., Eqs. 2.1.4 to 2.1.6). The mathematical form of linear filter 
G^(s) is given as follows: 


G 
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(s) 


G (s) 

V 



* VhJ 

+ — -I 


’ S + 
(S 


(3.1.1) 


(3.1.2) 
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► Velocity 


Figure 6 Modified Gaussian Turbulence Simulation 
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O(3)=0 


S + ~ A 

/3 


(S 4- v^/y 2 


where <|)q Is the white noise power spectrum. 

The linear filter, described above, is modified to include random 
variation of rms intensities. Random numbers generated by A are 
passed through a distribution modifier to generate rms intensities. 
Time histories are then generated by passing Gaussian white noise, <1)^, 
through the linear filter modified by the distribution modifier. 

I 

The patchy nature of atmospheric turbulence suggests that the 
turbulence field is composed of two components. One to represent 
variation of intensity within a patch and the other to represent 
variation of intensity with time (or from patch to patch) . The 
distribution modifier in this model, essentially, represents the 
variation of Intensity with time. The level of turbulence within 
each patch is controlled by the magnitude of the rms intensity, 

The distribution modifier Is the probability density function of 
the rms intensity. Analysis of several sets of atmospheric data 
characterized by various atmospheric conditions show that a truncated 
Gaussian distribution best fits the probability density of rms 
Intensity. (1) 


rms Distribution Modifier; 


1 r a - m 2 l 


where 
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P =■ probability density function 
“ rms intensity 

S = root mean square of rms intensity 
m = mean of rms intensity 
i = u, V, w gust components. 

Equation 3.1. A is completely described by the mean, m, and the root 
mean square, S, of the rms intensity. These variables have been 
derived from the data presented in Reference (1) characterized by 
terrain, altitude, and atmospheric stability. Table 1 represents the 
distribution modifier for two sets of atmospheric conditions. Through- 
out this report, the turbulence generated by these two distribution 
modifiers will be referred to as Model 2 and 3 (Model 1 is Gaussian 
turbulence simulation) . 

The scale lengths for these models are given by the Dryden form: 
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TABLE 1 



DISTRIBUTION MODIFIERS (rms INTENSITY) 


Mean Variance 


rms Distribution Modifier 
Model 2 


Altitude: 

250 ft C76.2 m) 

a 

u 

ft/sec 

(m/sec) 

3.1 

CO. 94) 

1.2 

CO. 37) 

Atmospheric 

Stability: 

Unstable 

a 

V 

ft/sec 

Cm/ sec) 

3.2 

CO. 97) 

1,2 

CO. 37) 

Terrain: 

Plains 

0 

w 

f t/sec 

(m/ sec) 

2.8 

CO. 85) 

0.9 

CO, 27) 


rms Distribution Modifier 


Model 3 


Altitude: 

750 ft C228.6 m) 

a 

u 

f t/sec 

Cm/sec) 

3.2 

CO, 97) 

0.8 

(0.24) 

Atmospheric 

Stability: 

Unstable 

a 

V 

f t/sec 

Cm/ sec) 

3.5 

(1,07) 

1.0 

(0,30) 

Terrain: 

Mountain 

a 

w 

ft/sec 

Cm/sec) 

4.1 

(1.25) 

0.9 

(0.27) 
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turbulence Intensity is, a , is given by 

w 

a n ^ ^ 

P(a ) = — exp (- i— ) (3.2.1) 

^ C2 

where is one-half the expected value of • 

Using Dryden spectrum models of real atmospheric turbulence, the 
value of C has been estimated in Reference (4) to be 2.3 ft/sec (0.70 
m/sec) . 

The rms intensity of the longitudinal, u, and the lateral, v, gust 
components are obtained from the relation: 


0 2 
u 




L 

u 


L 


V 


L 

w 


(3.2.2) 


The scale lengths are given by Equations 3.1.5 to 3.1.7. This will be 
referred to as Model 4. 

3.3 Variable Length and Intensity (VI.I) Turbulence Model 

The VLI turbulence model includes, in addition to the rms 
distribution modifier, a scale length modifier. A block diagram of 
this model is presented in Figure 7. In addition to controlling the 
patchiness of the turbulence field, the time variations of scale length 
achieves numerical compatibility with the real atmosphere and further 
randomizes the simulation. 

The scale length distribution modifier is derived from data 
collected in the LO-LO-CAT Program (1) for various combinations of 
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altitude, terrain, and atmospheric stability. Figures 8 and 9 show 
the fitted Gaussian distribution of scale length modifier for two sets 
of atmospheric conditions. The scale length distribution modifier is 
assumed to have the form 



where 

P = probability density function 

ttl 

= scale length of i component 
S = root mean square of scale length 
m • mean of scale length distribution 
i = u, V, w gust components 


(3.3.1) 


Table 2 presents the root mean square and mean of scale length 
distribution along with the rms distribution modifier for specific 
atmospheric conditions. The turbulence signal generated by these two 
atmospheric conditions will be referred to as Models 5 and 6. 



Probability (10~3) 
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0 Longitudinal Connponent^ 


O Lateral Component 
X Vertical Component 



S3Tnbols represent empirical data. 


Figure 9 Scale Length Distribution 
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DISTRIBUTION MODIFIERS (SCALE LENGTH) 
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Scale Length Modifier ^ ft to) 460(140.2) 126.6(38.6) 

Model 6 V 

L ft (m) 425 (129.5) 132.9 (40.5) 



CHAPTER IV 


THEORETICAL ANALYSIS OF PROPOSED MODELS 

In this section results obtained by statistical analysis of the 
gust velocity components for each of the six models will be discussed 
and compared with the properties of real atmospheric turbulence where 
possible. The statistical results have been obtained in the form of; 

1) mean and standard deviations 

2) normalized fourth and sixth moments 

3) probability density functions 

4) power spectral densities 

5) patchiness 

6) frequency of element of surprise. 

Table 3 tabulates the mean and standard deviation of gust 
components for each of the six models. It may be observed that the 
standard deviation varies from 2.6 to 5.2 ft/sec (0.79 to 1.58 m/sec) 
which is typical of low altitude clear air turbulence. 

Fourth and sixth moment characteristics are tabulated in Table 4. 
Within the limits of experimental error these characteristics for the 
VLI models are in fairly good agreement with the real atmospheric data 
obtained in Reference (5) . 

Since the cumulative probability and the probability density 
function essentially contains Identical information, only the 
probability density function will be analyzed. Figures 10 to 15 are 
plots of probability density functions for the simulated cases. In 
order to compare these with real atmospheric turbulence, a Gaussian 
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TABLE 3 


MEAN AND STANDARD DEVIATION OF GUST COMPONENTS 
CIO -min. sample) 


Gust Component 


Model 

No. 

Output 

Statistics 

u ft/sec 
(m/sec) 

v ft/sec 
(m/sec) 

w ft/sec 
(m/sec) 

Simulation 

Technique 

1 

Mean 

St, Deviation 

0.08 

3.97 

(0.02) 

(1.21) 

0.06 

3.90 

(0.02) 

(1.18) 

-0.03 

4.43 

(-0.009) 

(1.35) 

Gaussian 

9 

Mean 

0.83 

(0.25) 

-0,32 

(-0,09) 

-0.15 

(-0.04) 

Modified 


St. Deviation 

3,90 

(1.18) 

3.50 

(1.06) 

2.60 

(0.79 

Gaussian 


Mean 

0,88 

(0.27) 

-0.40 

(-0.12) 

0,06 

(0.02) 

Modified 


St. Deviation 

3.90 

(1.18) 

3.90 

(1.18) 

3.80 

(1.15) 

Gaussian 

/. 

Mean 

-0.36 

(-0.11) 

-0,16 

(-0.049) 

-0.22 

(-0.06) 

Rayleigh 

H 

St, Deviation 

5,19 

(1.58) 

4.84 

(1.47) 

4.48 

(1.36) 

Model 


Mean 

0.27 

(0.08) 

-0,36 

(-0.11) 

-0.20 

(-0.06) 

VLI 

J 

St, Deviation 

3.66 

(1.11) 

3,55 

(1.08) 

2.67 

(0.81) 

Model 

A 

Mean 

0.20 

(0.06) 

-0.10 

(-0.03) 

-0.33 

(-0.10) 

VLI 

VJ 

St. Deviation 

3.67 

(1.12) 

3.90 

(1.18) 

3.81 

(1.16) 

Model 



30 


TABLE 4 


NORMALIZED FOURTH AND SIXTH MOMENT DATA 


OF REAL AND SIMULATED TURBULENCE FIELDS 




(Over a 10-miii, sample) 


Model 

Normalized 

Gust 

Velocity 

Component 

Simulation 

No. 

Moment 

u 

V 

w 

Technique 

Real 

Fourth 

3.5 

3.5 

3.5 

Real atmospheric 

Atm. 

Sixth 

21,7 

21.7 

21.7 

turbulence data 


Fourth 

3,0 

3,0 

3.0 


1 

Sixth 

15.0 

15.0 

15.0 

Gaussian 

0 

Fourth 

5.9 

3.5 

3.2 

Modified 

Z. 

Sixth 

61.0 

22.3 

16.8 

Gaussian 


Fourth 

5.1 

3.2 

2.8 

Modified 

•J 

Sixth 

46.7 

18.9 

11.9 

Gaussian 

A 

Fourth 

3.7 

3.2 

3.3 

Rayleigh 

■t 

Sixth 

21.7 

18.1 

19.9 

Model 


Fourth 

3.5 

3.2 

3.5 

VLI 


Sixth 

20.8 

16.0 

21.8 

Model 

A 

Fourth 

3.1 

3.2 

3.9 

VLI 

O 

Sixth 

14.0 

16.1 

21.5 

Model 
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Figure 11 Probability Density of Simulated Field 



Probability % 
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Figure 12 Probability Density of Simulated Field 
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Figure 13 Probability Density of Simulated Field 






Model 6: VLI 


X Simulated Data 



3.58) 


Figure 15 Probability Density of Simulated Field 
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distribution is plotted on the same scale. It has been established (5) 
that real atmospheric turbulence exhibits a higher probability of both 
smaller and larger gust velocities compared to a Gaussian distribution. 
A careful study of the probability density of the simulated field 
reveals a higher probability of larger gust velocities compared to a 
Gaussian distribution, however the distributions, with the exception 
of Model 6, do not show higher probability of lower gust velocities. 

Power spectral densities of the simulated turbulence models are 
presented in Figures 16 to 21. The higher frequency components are 
compared with a line of slope -2 which is a characteristic of real 
atmospheric turbulence. The power spectrum in the entire frequency 
range within the limits of experimental error is in fairly good agree- 
ment with the assumed Dryden form (Equations 2,1,4 to 2.1,6), 

The patchiness of each of the models is plotted in Figures 22 to 
24, The derivative of vertical gust component is. plotted Illustrating 
a varying Intensity of patchiness. Model 6 presents patchy character- 
istics which closely match real atmospheric turbulence. 

Element of surprise Is tabulated in Table 5. At present there is 
no criterion available to either quantitatively measure this phenomenon 
or to establish a basis of comparison. In this report, "sudden jump" in 
the velocity field is used to describe element of surprise. 



Spectral Density (ft/sec)^/(rad/sec) 


100 


Model 1; Gaussian 



Frequency - CPS 


Figure 16 Power Spectral Density of the Simulated Field 
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Model 2: Modified Gaussian 



Frequency - CPS 

Figure 17 Power Spectral Density of Simulated Field 
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Frequency CPS 


Figure 20 Power Spectral Density of the Simulated Field 
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Model 6: VLI 



Figure 21 Power Spectral Density of Simulated Field 
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TABLE 5 


FREQUENCY OF ELEMENT OF SURPRISE 
OF SIMULATED FIELD * 

Frequency of Element of Surprise 


Model 

ja 

V 

w 

1 

0,03 

0.0 

0.0 

2 

0.07 

0.07 

0.0 

3 

0.03 

0.00 

0.0 

4 

0.0 

0.0 

0,23 

5 

0.03 

0.03 

0.0 

6 

0.27 

0.40 

0.0 


*For a 3.5 ft/sec Jump in velocity field). 



CHAPTER V 


TEST PROGRAM 

This chapter describes the flight simulator experiment, including 
the details of the aircraft simulated, the flight simulator, and the 
pilot performance task. 

5.1 Simulated Aircraft 

The aircraft simulated is the Canadian deHavilland DHC-6 Twin 
Otter. This particular aircraft is chosen as representative of light- 
wlng-loading STOL aircraft. In addition, there are pilots available 
with flying experience in the Twin Otter who can validate the simula- 
tion. 

Aerodynamic and dynamic stability parameters are listed in 
References (5) and (8). A summary is given in Table 6. 

5.2 Aircraft Simulator ^ 

The Visual Motion Simulator (VMS) at the NASA Langley Research 
Center, a synergistic motion-base simulator with the basic interior and 
Instrumentation of a jet transport cockpit (Figures 25 and 26), was 
employed in this study. A schematic diagram of the simulator, its con- 
trol system, and its data output capabilities is presented in Figure 
27 (8). A CDC-6600 digital computer, used exclusively to operate the 
real-time simulators, was programmed with the aircraft flight condi- 
tions, stability derivatives, six-degree-of-freedom differential equa- 
tions of motion, and a simulator washout routine. The program 

^This description has been adopted from Reference (8). 
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TABLE i3 

AIRCRAFT PARAMETERS (REFERENCE 8) 


w - 11500 lb (51152 N) I = 16900 slug-ft^ 

2 

Ug ” 256,67 ft/sec (78,2 m/sec) = 27600 slug-ft 

C„ ■= 0,045 I = 40600 slug-ft^ 

Og = -1.3° c = 6,5 ft (1,98 m) 

S •= 420 ft^ (39.0 m^) b = 65 ft (19.8 m) 

Longitudinal Stability Coefficients 

(22907 kg tn^) 
(37411 kg m^) 
(55031 kg tn^) 

C^ = 0.3818 

S. 

a 

= 0 


C = C - C^ 
X L D 

a < 

Cg => 0,045 

C 

m. 

a 

= -5.9 


C = -(L - 1 

z L 

a a 

C„ = 0.035 
M 

^L 

q 

= 5.504 


C = -C^ 

X. D. 

a a 

Cg = 5.7295 

s 

q 

= 0 


C = -CL 

z, L. 

a a 

Cg = 0.1432 
01 

C 

m 

q 

= -23.948 


C = -Cr, 

X D 

q q 

C = -1.9098 

c 

= -2 C_ 


C = -CL 

in 

a 

C =1.52 

1j « 

a 

X 

u 

c 

z 

u 

Lateral 

D 

= -' "l 

Derivatives 


z L 

q q 

C = -0.89 

c 

"p 

= -0.1 


C =0.5 

C = -0.12 

B 

P 

= -0.5488 


=0.13 

r 

C = 0.1215 

c 

n 

P 

Control 

= 0.006 
Derivatives 


C = -0.1855 
n 

r 

C = 0.39 

5 

= 0.0398 


C = -1.79 

"5 

C = 0.00348 

c 

"d 

= -0.1 


C ° = 0.45 

C ^ = 0.2055 
^<5 

a 

c 

"6 

a 

= -0.01 


e 


n a 
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Figure 25 The Visual Motion Simulator (VMS) at the NASA Langley Research Center (Ref. 8) 
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Figure 26 Interior Cabin of the Visual Motion Simulator (Ref. 8) 



Figure 27 Block Diagram of Motion-Base Simulator Apparatus (Ref. 8) 
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Integrates the equations of motion 32 times a second. These values are 
used by the simulator washout routine to determine the position of the 
simulator legs and the dynamic characteristics of the hydraulic 
actuators, the simulator is not capable of producing the magnitude and 
the duration of displacements, velocities, and acceleration of the real 
aircraft. It is the purpose of the washout routine to appropriately 
scale down the predicted motions of the real airplane to values that 
the simulator can produce without exceeding any of its design limita- 
tions. The washout routine also attempts to drive the simulator legs 
back to their neutral position following a disturbance from equilibrium 
in anticipation of a future disturbance, A detailed description of the 
physical dimensions and the performance specifications of the VMS may 
be found in References (9) and (10) . 

5,3 Pilot Task Performance 

Seven pilots, experienced in civil, military, and research flying 
were employed in the test program. Test runs, each of 8 to 10 minutes 
duration, for each of the six turbulence models were made in one pilot 
session. During separate sessions, some of the pilots repeated the six 
models in a random order. It was decided to have the pilots fly in a 
level flight constant altitude tracking task with no visual or "out-of- 
wlndow" cues in order not to introduce too many variables that might 
distract the pilots from their primary objective of trying to distin- 
guish differences between various turbulence models. After each run, 
the pilot was asked for his comment on the turbulence through use of a 
flight questionnaire (see Appendix B) Here the pilot was asked to 
estimate the turbulence intensity, realism, relative amplitude of 
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aircraft motions In each of the slx-degrees-of-freedom, patchiness, 
workload, task performance, and to give a Cooper-Harper handling 
quality rating (Figure 28) (11) for the airplane turbulence interaction. 
Additional questions explored the basis for the pilots' judgments. In 
addition, the pilots were also asked to estimate the altitude, terrain, 
and atmospheric stability in relation to his flying experience in 
turbulence. 

Several aircraft parameters, such as pitch, roll, yaw, and normal 
acceleration were recorded on strip charts for further analysis. A 
sample atrip chart is presented in Figures 29a and 29b. The rms 
intensities of the longitudinal, lateral, and vertical gust fields 
are presented in Table 7. The output intensity (rms) is the statisti- 
cal analysis of a ten-minute sample of the gust field. 



HANDLING QUALITIES RATING SCALE 
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Figure 28 Cooper-Harper Pilot Rating Scale 
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Figure 29a Typical Stip Chart Output of Flight Parameters 
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Figure 29b Typical Strip Chart Output of Flight Parameters 



lOTUT-OUTPUT INTENSITIES 
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CHAPTER VI 

RESULTS AND DISCUSSION OF SIMULATION 

Data obtained during the flight test program consisted of pilot 
opinion ratings and commentary relating to the simulated environment, 
aircraft handling qualities, and data relating to the physical 
environment to which the pilot was exposed. During each run, 
continuous strip chart recordings were made displaying time histories 
of various aircraft parameters for a later analysis. These include the 
three linear accelerations and three angular rates of aircraft in the 
body axes, elevator, aileron, and rudder deflections, throttle position, 
altitude, rate of climb, airspeed and aircraft heading. The pilot 
opinion ratings, obtained through a questionnaire, were in the follow- 
ing form; 

1) realism of turbulence 

2) correctness of relative amplitude of disturbances 

3) patchy characteristics 

4) frequency contents 

5) element of surprise 

6) atmospheric conditions 

7) handling quality ratings (Cooper-Harper) . 

These opinion ratings have been statistically analyzed and the 
results are presented in the following forms: 

a) Mean and Standard Deviation : of pilot opinion ratings 

for each of the turbulence models (Figures 30 to 37) . 

b) Correlation Matrix ; correlation among the various physical 
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Figure 30 Pilot Opinion Ratings 




Figure 31 Pilot Opinion Ratings 
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Figure 32 Pilot Opinion Ratings 
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Figure 33 Pilot Opinion Ratings 
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Ftgure 34 Pilot Opinion Ratings 
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Terrain 



Model Number 


A1 ti tude 



Model Number 


Figure 36 Atmospheric Condition 



stability 



Figure 37 Atmospheric Condition 
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characteristics of atmospheric turbulence is determined by using 

n _ 

I (x - x)(x - y) 

1 i=l 


xy n-1 


( 6 , 1 . 1 ) 


a a 
X y 


where 


xy 

x»y 


a fO 

X y 


correlation between x and y 
mean of x,y 

standard deviation of x,y 
number of observations. 


The correlation matrix is presented in Table 8. 

The following observations can be made from the statistical 
analysis of pilot opinion ratings: 

1) Figure 30 presents the pilot opinion ratings of handling 
quality and the realism of turbulence. It may be observed 
that at approximately the same rms intensity (see Table 7) 
of turbulence, the handling quality ratings transit from 
the satisfactory level, for a simple Gaussian model, to an 
unacceptable level for the more realistic and compositely 
structured VLI turbulence model. 

2) Figure 31 depicts the element of surprise and the patchiness 
ratings. The Gaussian model (Model 1) was found to be a 
little too continuous by almost all the pilots. On the 
other hand, Rayleigh model (Model 4) was rated "ahout right," 


as was Model 6. 
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3) Figure 32 presents the frequency content (low and high) 
ratings. The Gaussian model (1) was poorly rated whereas 
the mean ratings of Rayleigh and Vhl turbulence model were 
in the range of "about right." 

4) Figures 32 to 35 present amplitudes of disturbances as 
perceived by the pilots. The ratings show a progressive 
improvement as the pilots are exposed to more sophisticated 
models (see Models, 4, 5, and 6), 

5) Figures 36 and 37 present the atmospheric condition 
observations in the form of terrain, altitude, and 
atmospheric stability. The primary purpose of evaluating 
these was to determine how sensitive the pilots were to 
changes in atmospheric condition. Most pilots, when 
exposed to the six turbulence models, thought they were 
flying over level plains. On the altitude rating, the 
pilots flying the Gaussian model felt this turbulence 
was typical of altitude greater than 10,000 feet whereas 
they consistently rated the other models as typical low 
altitude turbulence. 

6) Table 8 presents the correlation matrix for various 
turbulence properties and aircraft handling qualities. 
Several Important observations can be made from this 
symmetric matrix. Realism of turbulence is highly 
correlated with patchiness (0.58), element of surprise 
C-0.63), and frequency content (0.52). This shows that 
in the opinion of pilots, the. realism of a turbulence 
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model Is closely linked to the physical properties of real 
atmosphere. In addition, the high correlation between 
handling qualities and realism CO, 74) indicates that the 
handling qualities are considerably worse for more realistic 
turbulence models. The low correlation between the patchi- 
ness characteristics and the intensity of turbulence (0.07) 
shows that the non-Gaussian patchiness characteristics cannot 
be Induced by simply chosing a higher level of intensity (rms) . 
On the other hand, patchiness is correlated to frequency 
contents (0,45) and the handling qualities (0.5), 



CHAPTER VII 


CONCLUSIONS 

This report has described several proposed turbulence models for 
producing artificial turbulence time histories which match the desired 
statistical properties of real atmosphere better than the presently- 
used simulation techniques. The use of these models gives Improved 
realism and accuracy in piloted simulator studies of handling 
qualities as affected by atmospheric turbulence. 

From the analytical study of the time histories generated by these 
models, and their comparison with real atmospheric turbulence, the 
following conclusions can be drawn: 

a) Turbulence simulated by the ’/LI gust models adequately 
matches the probability distribution (fourth and sixth 
normalized moments, probability density, and cumulative 
probability) of real atmospheric turbulence; and hence, 
presents an improved representation of atmospheric 
turbulence. 

b) Frequency content and the patchy characteristics of real 
turbulence can be closely matched. 

c) The proposed turbulence models (VLI) can flexibly 
accommodate changes in atmospheric conditions characterized 
by terrain, altitude, and atmospheric stability. This 
flexibility is not provided by any of the presently-used 
techniques . 
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d) The mechanization of the proposed models on a motion-base 
simulator is easy and inexpensive computationally because 
these models utilize only three linear filters. 

The time histories derived from turbulence models and the commonly- 
used Gaussian model were employed in a flight simulator experiment in 
order to determine the extent of pilot sensitivity to realism of 
various turbulence models and to evaluate the effect of turbulence on 
aircraft handling qualities. The principal conclusions drawn from the 
flight simulator study are; 

a) As expected from the analytical study, pilot opinion ratings 
show a considerable improvement in turbulence properties 
(realism, patchiness, frequency contents, etc.) over the 
most commonly-used Gaussian turbulence model, 

b) The correlation coefficient between the handling quality 

and the realism of turbulence is 0.74. This high correlation 
indicates that the handling qualities are considerably worse 
for more realistic turbulence models. 

c) From the flight test results of this program, it is apparent 
that the pilot's ability to handle the airplane in a turbulent 
environment not only depends on the rms intensity, but also 
the composition and the structure of turbulence. Pilots 
rated handling qualities in the satisfactory range while 
flying in a turbulence environment simulated by a simple 
Gaussian model; whereas the handling quality ratings degraded 
while flying in a turbulence environment simulated by the VLI 
turbulence model of approximately the same intensity. In 



fact, the handling quality ratings monotonically degrade 
as the pilots encountered more complex and realistic turbulence 
models. It may be concluded, therefore, that handling quality 
studies, using motion-base simulators, are critically affected 
by the suitable choice of a realistic turbulence model in 
addition to the appropriate rms intensities of turbulence. 

These tests were conducted in a simulated environment of a light 


general aviation STOL airplane. Caution should, therefore, be exercised 
in applying and extending the results to a general aircraft configuration. 
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APPENDIX A 


REVIEW OF BASIC DEFINITIONS C5) 

Statlonarlty ; A random process Is stationary if its statistical 
properties are not dependent on the time of their measurement. One 
could, for example, collect an infinite number of time histories, called 
an ensemble, which are representative of the process. If one takes an 
average across the ensemble, and if these averages are not a function 
of time, the process is stationary. 

Homogeneity : A random process is homogeneous if its statistical 

properties are Independent of position. 

Ergodiclty : In turbulence measurements it is impossible to obtain 

an ensemble from atmospheric measurements. Thus it is necessary to use 
time averages to get statistical information. If such a time average 
yields the same statistical properties as the ensemble average, the 
process is called ergodlc. 

Mean Value : The mean value of a random variable, u, of an ergodlc 

random process is given by 


u 


Lim 1 j-T 
X->«> 2T ■'-T 


(A.l) 


In practice the limit is not required and u can be approximated by 


u ^ Y /q u(t)dt. 


for T large. 


(A. 2) 


This approximate representation is especially useful for processes such 
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as turbulence. However, the time Interval T must be large enough so 
that the average approaches the asymptotic value one would obtain for 
a stationary process. 

Variance ; The variance of u is defined as 

ir ^-T " u)2jdt. (A. 3) 

As before in practical applications the variance can be approximated 
by 


*^u ~ Y -^0 “ ^^dt, for 


sufficiently large T. 


(A. 4) 


Standard Deviation (Root Mean Square ) : The standard deviation is 

defined as the square root of the variance. 

ttl 

Normalized Central Moment: The n normalized central moment, M 

— — _ — » n 

of a random process, u(t) , is 


M 


n 


Lira fT p u(t) - u 
X->« 2T -T ^ a 

u 


3 dt 


n = 1,2,3. . . (A. 5) 


which can be approximated by 


M 

n 



n = 1,2,3. . . (A. 6) 


Cumulative Probability Distribution ; The cumulative probability 
distribution of u(t) , is defined as the probability that u ^x. 
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Probability Density Distribution ; Probability density distribution 
of u(t) , P^Cx) Is defined as the probability that 

X < u ^ X + dx. 

Gaussian Probability Density Distribution ; If a random variable, 
u(t). Is Gaussian distributed its probability density is given by 

P^(x) = — exp C- Y (^Y ' ~ (A.* 7) 

a /liT °u 

u 

Rayleigh Distribution ; Another probability density of Interest 
is the Rayleigh distribution defined as follows: 

P(x) = — exp (- |— ) (A. 8) 

c2 2 ^2 

where Is one half the expected value of the random variable x or 

= Y E(x) = J /q xP^(x)dx (A. 9) 

Cross Correlation Function ; The cross correlation function of two 
random processes u(t) , w(t) is defined as 

*'■ (A. 10) 

correlations are the measures of the predictability of a signal at some 
future time (t + t) based on the knowledge of a signal at time t. 
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Autocorrelation Function ; The autocorrelation function is a 
special case of the cross correlation function defined above in which 
w(t) = u(t) , such that, 

\u^^^ = ^ ^-T 

Integral Scale Length : A statistical parameter of special 

importance in atmospheric turbulence is the integral scale length, 

u„ ^ 

L = — r R (x)dT, (A. 12) 

U 9 ■' -00 UU 

U 

where u^ is the reference steady state flight speed of the aircraft 
flying through turbulence. Scale length is an approximate measure of 
the distance an aircraft flies through turbulence. 

Cross Spectral Density : The cross spectral density of two random 

processes u(t) and w(t) is defined as the Fourier transform of their 
cross correlation 

$ (f) = / R (x)exp(- i 27 Tfx)dx , (A. 13) 

where f is frequency. 

Power Spectral Density ; The power spectral density, PSD, of a 
random process is the Fourier transform of Its autocorrelation function, 
or 
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({> (f) = / R (T)expCi2irfT)dT: . (A, 14) 

U ' IIU 

The PSD can be interpreted physically as the average contribution to 
the variable o^^from the frequency component f. Thus, 

a r (}) (f)df. (A. 15) 

White Noise ; White noise is a random process for which the PSD is 
a constant independent of frequency. That is, 

<J)^(f) = = constant, (A. 16) 



APPENDIX B 


FLIGHT QUESTIONNAIRE 

r I ir)ht Number Date 

Pilot: 

1. Turbulence Intensity: 

Li<jht Moderate Severe ___ Extreme 

2 . Realism of Turbulence: 

Very Good Good Fair Poor Very Poor 

3. Correctness of Relative Amplitude of Disturbances; 

Not Enough About Right Too Much No Comments 

Roll 


Pitch 

Yaw 

Heave ___ 

Side Force 

4. Patchy Characteristics (Variation of Intensity Bursts) 

Much Too Continuous A Little Too Continuous About Right 

A Little Too Patchy No Comments 

5. Frequency Contents of Turbulence; 

Not Enough About Right Too Much No Comments 

Low FRQ: 

High FRQ; 

6. Element of Surprise in the Simulated Turbulence Field; 

a. Quite Often Sometimes Never 

b. Real ism of 6a: 

Very C5ood Good Fair Poor Very Poor 


81 



82 


7. Atmospheric Conditions: 

a. Altitude: 0 -1,000 Ft I ,000 - 1 0,000 Ft 

Over 10,000 Ft Unable to Judge 

b. Atmospheric StabI I I ty: Stable Unstable 

Neutral Unable to Judge 

c. Terrain: Mountains Plains Unable to Judge 

8. Pilot Estimate of the Work Load: 

Very Easy Easy Average Difficult Very Difficult 

9. Pilot Estimate of Task Performance: (Integral Squared Error for ILS 

I raoki nq Task) 

Very (Tood Good Average Poor Very Poor 

ID. Realism of This Model Compared to Previously Flown Model: 

Very Good Good About the Same Poor Very Poor 

11. Did You Observe a Repetitive Pattern in the Turbulence Field?' 

Yes No 

12. Cooper-Harper Rating: 

13. Additional Comments About Realism of Turbulence and Aircraft. Simulation: 
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APPENDIX B 
(Cont.) 

PILOT EXPERIENCE 


I . Nonio Date 

7 . Wh<i1 Fypo of flying Experience Have You Had? 

Military Civil 

Main Types of Aircrafts Flown; 


4. fofal Number of Hours Flown: 

b. Hours of Instrument Flying: 

(). Hours in Simulators; 

/. Hours in VMS; _____________ 

H. Hours in Twin Otter: 

9. a. Estimate ttie % of Time Flown in Turbulence; 

b. Of This Time What ^ Was Flown in 

Light Turbulence Moderate Turbulence Severe Turbulence Extreme Turbulence 


lO. What Characteristic of Turbulence Interferes Most with Your Ability to 
Control the Aircraft? 


It. Describe the Most Critical Case of Turbulence Encountered During Your 
flying Experience: 

a. Day Night 

b. Terrain: Altitude: 


-■ PAGE IS 

■ ytT'R QUA.Ln^ 



84 


c. Atmospheric Stability: 

Stable Neutral Unstable Unable to Judge 

(I . What Was the Task You Were Attempting Before Turbulence Was Encountered: 

(e.g. I LS Approach, Cruise, etc.) 

e. Any Additional Comments: 


ORIGINAL PAGE IS 
OF POOR QUALITY 





